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Resume —Cette etude propose une technique de caracterisation thermoaeraulique de mini-regenerateurs; organes clefs des mini- 
refroidisseurs croygeniques Stirling ou a tube a gaz pulse employes en refroidissement de capteurs de rayonnement, de detecteurs 
infrarouges ou de composants electroniques. Nous utilisons une adaptation de la methode en regime transitoire, dite « single-blow », 
qui consiste a perturber subitement I’equilibre thermique du regenerateur en le faisant traverser par un fluide chaud de temperature 
etablie a son entree. La modelisation thermique repose sur un modele direct analytique convectif a trois equations pernant en compte, 
outre les classiques echanges de chaleur internes matrice/fluide, les pertes thermiques au travers de I’enveloppe du regenerateur et 
la convection naturelle externe conditionnant les performances thermiques dans les mini-systemes. L’estimation de trois coefficients 
d'echange de chaleur est obtenue en resolvant un probleme inverse de traitement numerique de donnees bruitees en temperatures 
issues de micro-thermocouples. Les resultats dynamiques et thermiques sont presentes pour un regenerateur tubulaire de petite 
taille rempli de billes en acier, dans lequel s’ecoule de I'air a des valeurs du nombre de Reynolds hydraulique comprises entre 200 et 
1 350. © 2001 Editions scientifiques et medicales Elsevier SAS 

mini-refroidisseur / mini-regenerateur / methode transitoire / probleme inverse / coefficient d’echange de chaleur / 
micro-thermocouple 

Abstract —Characterization technique of a thermal mini-regenerator for Stirling or pulse tube mini-refrigerator. This study 
proposes a technique of dynamic and thermal characterization of mini-regenerators; the most important piece of Stirling and 
Pulse Tube cryogenic refrigerators used for cooling of radiation sensors, infrared detectors or electronic components. We use an 
adaptation of the transient method known as “the single-blow technique” which consists in disrupting suddenly the regenerator 
thermal equilibrium by a heated fluid flowing through it with an established level of inlet temperature. The thermal modeling is based 
on a three-equation analytical convective model taking into account, in addition to the traditional internal heat transfer between 
matrix and fluid, heat losses through the regenerator wall and external natural convection which govern heat transfer performance 
in mini-systems. The estimation of three heat transfer coefficients is obtained by resolving an ill-posed inverse problem with noisy 
temperature measurements from micro-thermocouples. The dynamic and thermal results concern a regenerator matrix made of steel 
spheres inside a small circular tube with air flowing at hydraulic Reynolds numbers ranging from 200 to 1 350. © 2001 Editions 
scientifiques et medicales Elsevier SAS 

mini-refrigerator / mini-regenerator / transient method / inverse problem / heat transfer coefficient / micro-thermocouple 


Nomenclature 

A conductance hydraulique . . . . kg-s - 1 -p a ~l 

a coefficient de la loi de frottement 

b coefficient de la loi de frottement 
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B, S 0 , Si, 

S2, S e coefficients du modele 

Bi nombre de Biot 

C[ coefficient de frottement de 

Darcy-Weisbach modifie 

c capacite calorifique. J-kg _ 1 -K _1 

diametre hydraulique. m 

D\ diametre interieur du regenerateur m 

d diametre de til d'une grille ... m 

e epaisseur du tube. m 
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debit massique. kg-s 1 

coefficients de convection 

surfacique. 

permeabilite intrinseque de la 

matrice . nr 

nombre de Knudsen 
longueurs, respectivement du fil et 

du regenerateur. m 

dimension caracteristique ... m 

fibre parcours moyen des 

molecules de gaz. m 

masse molaire du gaz. kg-mol _1 

variable de Laplace 
adimensionnelle 

pression absolue locale. Pa 

rayon polaire. m 

constante des gaz parfaits . . . J-K _ 1 mol _1 

temps caracteristique du 

regenerateur. s 

nombre de Reynolds hydraulique 
nombre de Reynolds de particule 


section totale. 

temps. 

temperature absolue locale . . . 
temperature absolue moyenne . 

vitesse de Darcy. 

vitesse interstitielle. 

coordonnee d’espace. 


m 


2 


s 

K 

K 

m-s -1 

m-s -1 


m 


Symboles grecs 

/x viscosite dynamique. kg-ni - -s —1 

p masse volumique. kg-m -3 

&7Z ■> 

02, o e surface specifique. m _1 

e porosite 

kq coefficient de forme 

X conductivity thermique. 

y constante isentropique d'un gaz 

r temps adimensionnel 

9 temperature reduite 

X coordonnee d’espace reduite 

Indices et exposant 


a ambiance 

e exterieur ou entree 

i interieur 

g gaz 

s solide 

t tube enveloppe 

fib surface fibre de passage du fluide 

y indicateur simplificateur d’ecriture 

imp impulsionnelle 


1. INTRODUCTION 


Les regenerateurs ont ete developpes en 1816 par Ro¬ 
bert Stirling pour la mise au point de son moteur a air 
chaud. 11s represented une importante classe d’echan- 
geurs de chaleur dans lesquels les transferts thermiques 
sont obtenus par F intermediate d’un echange periodique 
entre un gaz et une matrice solide. En fait, le role du 
regenerateur pourrait etre ramene a celui d’une eponge 
thermique : prendre et restituer alternativement de la cha¬ 
leur au gaz le traversant pour le refroidir de la tempera¬ 
ture ambiante a celle de la partie froide. Ainsi, lorsque 
le gaz va du compresseur vers l’echangeur de la source 
froide, il est «thermalise» par le regenerateur. C’est un 
echangeur absolument crucial a la realisation des cycles 
thermodynamiques dans les refrigerateurs cryogeniques 
de types Gifford-MacMahon, Stirling et a tubes a gaz 
pulse (TGP). L’optimisation des performances frigori- 
fiques dans le cas de microrefroidisseurs qui doivent al- 
lier haute efficacite et petite taille passe par les connais- 
sances precises des transferts thermiques internes et de la 
perte de charge. Dans ces systemes, le regenerateur est le 
plus souvent constitue d’un material! poreux sous forme 
d’un assemblage de petites grilles, de tissus metalliques 
fins ou de microbilles. Au cours des cycles successifs, 
les variations de la temperature de la matrice solide doi¬ 
vent etre reduites au maximum en utilisant un empile- 
ment de longueur importante a haute capacite thermique, 
dans lequel la conduction axiale doit etre la plus faible 
possible en raison de la position du regenerateur entre 
l’organe de compression et la source froide du refrige- 
rateur. Le volume de gaz contenu dans le regenerateur 
represente un volume mort supplementaire pour le mini- 
refroidisseur et penalise done ses performances. La poro¬ 
site et la longueur de la matrice deviennent alors des pa- 
rametres importants pour lesquels il faut adopter un com- 
promis en raison de l’impossibilite de concilier efficacite 
frigorifique elevee et perte de charge faible. A ce titre, 
le remplissage par billes amene des porosites inferieures 
au remplissage par grilles (porosite e 40 % au lieu de 
60 %) et par contre coup des pertes de charge plus impor- 
tantes. Les propos precedents montrent bien la difficulte 
de concilier les differentes fonctions du regenerateur. Le 
meilleur compromis en matiere de caracterisation com¬ 
plete des performances d’un regenerateur donne ne peut 
done sortir que d’une analyse globale prenant en compte 
a la fois la dynamique des fluides et les echanges ther¬ 
miques. La difficulte est accrue en raison du fonctionne- 
ment particulier de ces machines : debit de gaz, pressions 
et temperatures internes variables dans le temps, thermo- 
dependance des proprietes physiques du fluide et des ma- 
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teriaux en raison de 1’ important gradient de temperature 
le long du regenerateur... 

Le dimensionnement actuel des regenerateurs indus- 
triels repose sur les travaux anciens de Nusselt, les 
abaques d’efficacite de Hausen ou le concept de nombre 
d’unite de transfert (NUT) de Shah [1], Malheureuse- 
ment, ces methodes sont mal adaptees aux mini-regene- 
rateurs en raison de leur imprecision eu egard aux condi¬ 
tions de fonctionnement tres particulieres des cryorefroi- 
disseurs. 

La litterature fait mention de trois methodes de carac¬ 
terisation des regenerateurs : 

• la methode calorimetrique; 

• la methode statique dite « single-blow » en anglais 
[2-4]; 

• la methode dynamique ou alternee [5, 6], 

Pour la presente etude, nous avons opte pour une adap¬ 
tation de la seconde en soumettant brusquement le rege¬ 
nerateur de petite taille, initialement isotherme et sans 
ecoulement, a un pseudo-echelon de temperature d’en- 
tree du gaz. Ainsi, F etude des performances d’un petit 
regenerateur tubulaire rempli de billes d’acier aleatoire- 
ment disposees et enclin aux pertes thermiques est abor- 
dee en utilisant Pair comme fluide de travail. Les resultats 
globaux obtenus portent tant sur le plan de la recherche 
de la courbe de frottement intrinseque a Forgane etudie 
que sur la description du transfert de chaleur a partir d’un 
modele convectif a trois equations (modele direct). La re¬ 
solution d’un probleme inverse de traitement numerique 
de donnees bruitees en temperatures issues de microther¬ 
mocouples permet de remonter, via Falgorithme de mi¬ 
nimisation de Simplex, aux valeurs de trois coefficients 
d’echange de chaleur globaux caracteristiques du rege¬ 
nerateur etudie sous les conditions experimentales impo- 
sees. La rapidite est Favantage principal de la technique 
de caracterisation developpee puisqu’elle repose sur un 
modele de connaissance analytique et demande peu de 
donnees experimentales a inverser pour obtenir des re¬ 
sultats pratiques. En absolu, seuls les releves temporels 
des temperatures du gaz d’entree et de sortie suffisent 
pour estimer les trois coefficients d’echange du modele. 
En particulier, les distributions spatiales des temperatures 
dans les billes et dans le tube a differents instants ne sont 
pas necessaires. De plus, la methode iterative d’estima- 
tion parametrique devrait apporter une meilleure preci¬ 
sion que les schemas classiques de dimensionnement em¬ 
ployes jusqu’a maintenant consistant a traiter les points 
experimentaux avec la methode des moindres carres or- 
dinaires ou a utiliser la methode de la tangente (maximum 
slope-technique) [2-4], 


2. MODELISATION DU PROBLEME 

2.1. Domaine de validite des equations 
classiques 

Le systeme complet d’equations pour un regenerateur 
miniature contenant un milieu poreux et traverse par un 
flux gazeux est normalement decrit a partir des relations 
generates locales de la mecanique des fluides et de la 
thermique. II comporte les trois equations de quantite 
de mouvement pour le fluide, F equation de continuite et 
les trois equations de l’energie pour le gaz, la matrice 
et le tube enveloppe. Un aspect fondamental qu’il faut 
envisager pour les systemes a echelle reduite traverses 
par un ecoulement gazeux est la validite des formulations 
classiques. Pour les ecoulements de fluide, on sait qu’il 
y a passage du regime continu de Poiseuille au regime 
moleculaire de Knudsen des que le libre parcours moyen 
des molecules de gaz L p devient important face aux 
dimensions caracteristiques L c du systeme. Le critere est 
le nombre de Knudsen : 

Kn = ^ « 1 (1) 

L p est calcule par la theorie cinetique des gaz via la 
relation suivante : 

T InM 

L p = —3L g (v — 1),/- (2) 

p P s r V 8 TR 

Dans la situation envisagee, le libre parcours moyen est 
de l’ordre de 10 -7 m; alors que la dimension des pores de 
la matrice avoisine 10 -4 m. On peut done conclure que 
les modeles habituels sont encore valables (Kn ~ 10 -3 ). 

2.2. Hypotheses utilisees 

La complexite des conditions experimentales rencon- 
trees ne permet pas d’effectuer une resolution analytique 
du probleme si des simplifications ne sont pas introduites 
par la suite. Les hypotheses suivantes sont retenues : 

• la masse volumique du fluide est consideree constante : 
cela suppose que le gradient de pression le long du 
regenerateur doit etre faible; 

• Fecoulement monodimensionnel est axial (absence de 
tortuosite): cette hypothese est verifiee si le rapport entre 
de diametre interieur du regenerateur I), et le diametre de 
particule d p est tel que Di/d p > 10; 

• la porosite de la matrice est uniforme (on retrouve la 
condition precedente D\/d p > 10); 
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• le fluide est assimile a un gaz parfait; 

• les transferts thermiques avec pertes sont caracterises 
par trois coefficients globaux de convection: gaz/matrice, 
gaz/paroi et paroi/milieu exterieur; 

• la conduction thermique axiale dans la matrice et le 
tube est negligeable : il faut pour cela que le gradient de 
temperature soit negligeable au sein de la matrice solide 
et du tube en satisfaisant la condition suivante sur le 
nombre de Biot: 

h y dry 

lily = ~A~ « 1 (y = S, t) 

A Ay 

• les variations des proprietes physiques sont negli- 
geables dans la plage de temperature etudiee; 

• l’effet de la detente Joule-Thomson est neglige dans 
nos conditions operatoires. 


2.3. Caracterisation dynamique 
du regenerateur 


De facon analogue aux tuyaux cylindriques, on carac- 
terise globalement la chute de pression motrice A p d’un 
ecoulement permanent d’un fluide visqueux en regime la- 
minaire ou turbulent a travers un tube rempli d’un mate- 
riau poreux a l’aide du coefficient de frottement modifie 
C| de Darcy-Weisbach : 


C't = 


2ApDh 

PgU 2 L 


(3) 


Dans le cadre de F etude des milieux poreux composes de 
microcanaux, des correlations de calcul sont proposees 
sous la forme suivante; a et b dependant du type de 
matrice : 

, a 

C' f =— + b (4) 

Re h 

Parmi la grande variete des definitions du nombre de Rey¬ 
nolds, la plus realiste decoule de F analyse de F ecoule¬ 
ment en conduite avec le diametre hydraulique I)\ } et la 
vitesse interstitielle du fluide u dans les pores de la ma¬ 
trice comme grandeurs de reference : 


Re h = 


PgU D}) 

Mg 


Cm D], 

/XgiSlib 


(5) 


Si le regenerateur est constitue de billes ou de grilles 
metalliques, le diametre hydraulique de matrice sera 
respectivement calcule par les relations suivantes : 


et la porosite e pour un grillage carre vaut: 

,=i- " d ji+(—X 

4 (l + d)]j \l + dj 

Le debit massique G m a pour expression : 

Giji PgM >Sjib Pgll£ S —— pg U S 

On definit alors la conductance hydraulique A : 

A G m 2spgD^S 

A p RehCf^igL 


(7) 

( 8 ) 

(9) 


On remarque que la relation entre la perte de charge, ou 
le gradient de pression, et le debit massique s’ecarte rapi- 
dement d’une simple loi lineaire puisque la conductance 
hydraulique A depend du debit massique par l’interme- 
diaire du nombre de Reynolds. 

Pour un empilement non consolide de particules iden- 
tiques simples, dans lequel transite un ecoulement lami- 
naire incompressible de Darcy, la loi de Kozeny-Carman 
permet de calculer la permeabilite K par : 


K = 


d 2 e 3 


36at 0 (1 - e) 2 


pour Re h <10 


_Dh 

Sdn 


GO) 


kq est un facteur de forme (3,6 < kq < 5). 11 vaut 4,8 pour 
un empilement de spheres de diametre c/ p . 

Ainsi, on montre que les coefficients a et b necessaires 
au calcul du coefficient de frottement modifie (equa¬ 
tion (5)) ont pour expression : 


2 _ 

K 


et b — 0 


( 11 ) 


Pour un ecoulement turbulent ou la degradation d’energie 
est plus importante, il faut utiliser la loi d’Ergun : 


a — 300 


(1 - g) 2 

e 2 d 2 


(1 -s) D h 

b — 3,5--- (12) 


Quelque soit le regime d’ecoulement isotherme, la 
determination experimentale du coefficient de frottement 
modifie C ( ' ne necessite que des mesures de pression et la 
connaissance du debit massique. 


2.4. Modele de transferts thermiques 


2 

D h =- 


edr i 


3 (1 — e) 


ou Dh = 


sd 


(1 - e) 


(6) 


Il s’agit de modeliser le comportement thermique en 
regime transitoire de regenerateurs de petite taille soumis 


458 



Technique de caracterisation d’un mini-regenerateur thermique 


a des pertes de chaleur avec l’exterieur non negligeables 
devant les echanges gaz-matrice internes au systeme. Par 
consequent, nous introduisons trois coefficients globaux 
de convection : le coefficient h e de pertes externes de 
Fenveloppe, le coefficient d’echange interne gaz/paroi /z t 
et le coefficient d’echange /? s entre le gaz et la matrice. 
Conformement aux hypotheses donnees precedemment, 
les trois equations locales de la chaleur s’ecrivent: 


Py c y 


8Ty 

dr 


Uy 3 Ty 
s dx 


= A.* 


d z T v 1 dTy 


dr 2 r dr 

avec U s = U\ — 0 


(y — g> s, t) (13a) 
(13b) 


Les conditions aux limites et initiale du probleme 
sont: 

• en r — 0 : les champs de temperature du gaz et de la 
matrice sont axisymetriques; 


97g 

dr 


— 0 

r=0 


et 


97; 

dr 


=0 

r=0 


(14) 


• en r = r; : on considere l’identite du flux thermique 
transfere; 


9?g 

’ dr 


= X t 


dT t 

dr 


et A s 


97s 

dr 


= 0 (15) 


• en r — r e : les pertes thermiques externes par convec¬ 
tion naturelle sont prises en compte; 


97) 

—- 

dr 


= h e (T a - T t \ r=re ) 


(16) 


La condition initiale d’uniformite des temperatures au 
sein de la totalite du regenerateur est imposee : 


Ty (0, x,0)—T a (y = g,s,t) (17) 


Apres moyennage spatiale des equations locales de la 
chaleur selon le rayon, le probleme thermique est ecrit 
sous forme adimensionnelle en introduisant les grandeurs 
suivantes : 


X = 

r 


A 

L 

t 

~Rc 


Ty(X, t) - T & 
Ty e ~ T-d 


(y = g. s, t) 


L p s c s ( 1 - e) + ptc t (2 + e/rj)e/n 
U PgCg 


(18) 

(19) 

( 20 ) 
( 21 ) 


B = R c 


h sCTp 


B i = 


(1 — e)p s c s 
L so\h[ 

U pgCg 


O e hr 

= — R c 

p t C t 


Bo 

B 2 


h s o p L 
PgCg U 


p t c t 


R c 


( 22 ) 


Les coefficients constants B , Bo, B\, B 2 et B e sont des 
expressions analytiques des parametres physiques et geo- 
metriques qui caracterisent les differents echanges ther¬ 
miques [7], R c represente un temps caracteristique qui 
situe Fabscisse du point de tangence a la courbe de la 
temperature de sortie du fluide au voisinage du point d’in¬ 
flexion. 


Dans Fespace de Laplace, la solution analytique du 
probleme thermique instationnaire prend la forme sui- 
vante : 


C p 0lb*) 

= exp ( ~ P 


Le 


Bo 


R C U p* + B 
P* + B e Bi 

p* p* + B e + B 2 


B9, 


imp 


r p (x,B*) = ^ 
fTW*) = 


B 


B 2 0, 


imp 


B e + B 2 


(23) 

(24) 

(25) 


Les reponses impulsionnelles 7’ imp (/,r) dans Fespace 
temporel sont obtenues par inversion numerique des 
transformees de Laplace e nnp (x , p*) (algorithme d’in- 
version numerique de Stehfest [8]). Les reponses en tem¬ 
perature pour une temperature d’entree du fluide quel- 
conque Tge (0, r) sont obtenues par le calcul numerique 
du produit de convolution suivant: 


7’(x> t) = T ge (0, r) • 7’ imp (/, r) (26) 


2.5. Technique de caracterisation 
thermique du regenerateur 

La technique de caracterisation thermique du re¬ 
generateur consiste a determiner les trois coefficients 
d’echange a partir des seuls releves temporels experi- 
mentaux des temperatures d’entree (echelon de tempe¬ 
rature) et de sortie du gaz. II s’agit de resoudre un pro¬ 
bleme inverse d’optimisation parametrique visant a faire 
correspondre un modele analytique avec des resultats ex- 
perimentaux. Ainsi, une procedure iterative minimise la 
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Figure 1 . Schema de I'installation experimental. 
Figure 1. Sketch of the test apparatus. 


somme des ecarts quadratiques entre les valeurs mesu- 
rees des temperatures de sortie du fluide et celles recalcu- 
lees par le modele pour fournir des estimateurs des coeffi¬ 
cients d’echange. L’algorithme d’optimisation utilise (al- 
gorithme de Simplex, Matlab) est particulierement adapte 
au traitement numerique de donnees experimentales tres 
bruitees, telles que celles delivrees par des sondes ther- 
moelectriques de petite taille. Dans une etude prelimi- 
naire [7], la presente technique a ete validee sur des don¬ 
nees simulees par le modele et artificiellement bruitees. 


3. DISPOSITIF EXPERIMENTAL 

L’element principal de I’installation experimentale, 
representee a la figure 1, est le regenerateur thermique 
tubulaire a billes. II se presente sous la forme d’un tuyau 
cylindrique isolant en polychlorure de vinyle (PVC), de 
diametre interieur D\ de 30 mm, de longueur L egale a 
94 mm et d’epaisseur e reduite a 5 mm pour minimiser la 
capacite thermique. L’isolation du tube enveloppe est vo- 


lontairement negligee puisque le modele prend en compte 
les pertes thermiques externes. Le remplissage est assure 
par des billes en acier, de diametre d p egal a 2,9 mm a 
deux pour cent pres. L’empilement aleatoire determine 
une matrice solide, de porosite s egale a 0,47 mesuree 
avec une incertitude de deux pour cent, assimilable a un 
milieu poreux non consolide. Les billes sont maintenues 
entre deux grilles metalliques, generant de faibles pertes 
de charge inferieures a 40 Pa, implantees respectivement 
dans le divergent d’entree et dans le convergent de sortie 
du regenerateur. Toutes les caracteristiques geometriques 
et physiques importantes du regenerateur et de l’empile- 
ment sont resumees dans le tableau I. 

Le principe de la methode de caracterisation du re¬ 
generateur a etudier consiste a determiner le coefficient 
de frottement modifie et les trois coefficients de convec¬ 
tion en reponse a un echelon de la temperature d’entree 
pour un debit de fluide constant. Le gaz utilise est Fair, 
dont les proprietes thermophysiques sont considerees va¬ 
riables avec la temperature selon les expressions du ta¬ 
bleau II. 
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TABLEAU I 

Proprieties geometriques et physiques (a 29B K) du regenerateur. 



P 

c 

>, 

dp 

£ 

L 

e 

A 

D[/dp 


kg-rrU 3 

J-kg-'-KT 1 

W-rru'-K -1 

m 

- 

m 

m 

m 

- 

Spheres (acier) 

7 660 

481 

52 

0,0029 

0,47 

- 

- 

- 

10,4 

Tube (PVC) 

1450 

1000 

0,16 

- 

- 

0,094 

0,005 

0,03 

- 


TABLEAU II 

Variation des proprietes thermophysiques de I’air en fonction de la temperature. 


Masse volumique 

Viscosite dynamique 

Conductivite thermique 

Capacite calorifique 

kg-rrU 3 

kg-rru'-s -1 

W-nW'-K- 1 

J-kg-'-KT 1 

p(p,T) = pM/RT 

p{T) = 1,875-10“ 5 (T/293) 0,6386 

\{T) = 0.026(7/293) 0,7231 

c g = 975,3+ 0.0368(7-273) 

+ 2-10^ 4 (7 — 273) 2 


Deux regenerateurs identiques sont montes en pa¬ 
rallel avec des parcours aerauliques absolument sem- 
blables de maniere a respecter la meme perte de charge 
dans chacune des branches. Un jeu de quatre electro- 
vannes pilotees par un microprocesseur et disposees en 
amont et en aval des deux echangeurs permet d’orienter 
le gaz vers le circuit requis de maniere strictement syn¬ 
chronise. Cet artifice permet de maintenir le debit de gaz 
desire et de «thermaliser »le gaz et le circuit avant de de- 
clencher F echelon de temperature d’entree du gaz, tout 
en conservant le regenerateur a tester a temperature ini- 
tiale uniforme (temperature ambiante). La mise en tem¬ 
perature du gaz pour realiser F echelon est assuree par un 
rechauffeur sous forme d’une resistance electrique enfer- 
mee dans un conduit dans lequel circule Fair. Une vanne 
pointeau permet d’imposer differents debits dans le rege¬ 
nerateur sous une pression constante d’entree ajustee au 
moyen d’un detendeur de precision. L’installation fonc- 
tionne en circuit ouvert. 

Les mesures instantanees de temperatures de gaz en 
regime transitoire imposent le choix de capteurs de pe- 
tites dimensions a faible inertie thermique pour suivre 
correctement les variations de temperatures. Pour ces rai¬ 
sons, notre choix s’est porte sur des microthermocouples 
Chromel/Alumel (type K) de diametres egaux a 12,7 pm 
a F entree et a 25 pm a la sortie, pour lesquels les valeurs 
des constantes de temps convectives ne depassent pas 
50 ms dans nos conditions operatoires. La technique de¬ 
licate de fabrication consiste a effectuer sous binoculaire 
un soudage des fils par decharge capacitive. Les deux 
microthermocouples sont implantes au cceur de l’ecou- 
lement a Fentree du regenerateur pour le 12,7 pm (pTl) 
et a sa sortie pour le 25 pm (pT2) dans le divergent et le 
convergent respectivement. Ces sondes sont etalonnees 


en regime thermique permanent dans des conditions ex- 
perimentales identiques a celles rencontrees lors des es- 
sais. Pour cela, on utilise un capteur de temperature de 
type semi-conducteur AD22103 1 , de sensibilite mesu- 
ree egale a 26,5 mV-K -1 et de grande fidelite dans la 
plage de temperature 273-373 K, monte successivement 
a F entree et a la sortie du regenerateur. 


4. RESULTATS EXPERIMENTAL^ 

4.1. Caracterisation hydrodynamique 
du regenerateur 

Nous avons effectue les essais de caracterisation hy¬ 
drodynamique du mini-regenerateur avec des debits al¬ 
lant de 1,3 a 8 Nnr’-h -1 sous une pression constante 
d’entree egale a 2,12 bars mesuree avec un manometre 
Bourdon (classe 0,5). L’incertitude sur la masse volu- 
mique est inferieure a un pour cent. Le comportement 
du regenerateur est etudie pour trois temperatures d’en¬ 
tree du gaz 7ge en regime «thermalise» (298 K, 313 K 
et 338 K). Les debits sont mesures avec une incertitude 
maximale de trois pour cent pour la pression de travail 
choisie. La mesure de la chute de pression au moyen 
d’un manometre a pression differentielle est obtenue avec 
une incertitude maximale de six pour cent pour les plus 
faibles debits. La contribution des grilles de maintien de 
l’empilement de billes n’excede pas trois pour cent de la 
chute de pression totale. Les proprietes thermophysiques 


1 Commercialise par la Societe One Technology Way, P.O. Box 9106, 
Norwood, MA 02062-9106, USA. 
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Figure 2. Evolution du coefficient de frottement modifie C' f en 
fonction du nombre de Reynolds pour une pression constante 
d'entree du gaz egale a 2,1 2 bars. 

Figure 2. Evolution of the modified friction factor versus the 
Reynolds number for a constant entrance gas pressure of 2.1 2 
bars. 


ont ete calculees aux differentes temperatures a l’aide des 
relations du tableau II. 

L’ensemble des resultats est presente a \& figure 2 sous 
forme de F evolution du coefficient de frottement modi¬ 
fie C{ en fonction du nombre de Reynolds Re h de plage 
de valeurs comprises entre 150 et 1 100. Dans nos condi¬ 
tions operatoires, les incertitudes relatives maximales sur 
les valeurs du nombre de Reynolds et du coefficient C' { 
sont estimees respectivement a 12 % et a 20 %. 

On observe que Fensemble des points experimentaux 
decrit une evolution tres voisine de celle qui se rapporte 
a la relation (12) (cas turbulent). La valeur du coefficient 
est alors comprise entre 2,01 et 4,33. 

Toutefois, nos mesures surestiment les resultats cites 
pour les faibles valeurs du nombre de Reynolds (jusqu’a 
+37,1 %). A Finverse, une sous-estimation du coefficient 
est relevee pour les fortes valeurs du nombre de Reynolds 
(jusqu’a —14,3 %). Un passage preferentiel du gaz a la 
paroi diminuant le coefficient de frottement n’est pas a 
exclure meme si le rapport DJdp est superieur a 10. La 
non uniformite en taille des particules solides et la nature 
compressible du fluide (considere comme isovolume) 
peuvent egalement expliquer ces ecarts. Notons aussi 
que dans nos conditions operatoires, la perte de charge 
peut engendrer une diminution maximale de la masse 
volumique de 3 %. 

Pour les ecoulements suffisamment turbulents, une 
faible dependance du coefficient de frottement aux va¬ 
riations du nombre de Reynolds est obtenue ici; il s’agit 
d’une tendance relevee par tous les auteurs de la littera- 
ture. 

On observe egalement une diminution significative du 
coefficient de frottement avec la temperature a nombre 
de Reynolds constant. Par exemple, si l’on compare 


les courbes de frottement relatives a T ge — 298 K et a 
7ge = 338 K, la diminution atteint onze pour cent pour 
Re h = 350. 

Au vue de Failure de la courbe relative a la tempera¬ 
ture 7g e = 298 K, une correlation de frottement a ete eta- 
blie avec un ecart relatif maximal de deux et demi pour 
cent par rapport aux mesures correspondantes : 


, 406,6 

C = --—F 1,76 

f Re h 


(27) 


Cette correlation presente des ecarts relatifs maximum 
de quatorze pour cent et trente trois pour cent avec la 
correlation de Ergun pour les grands et les petits nombres 
de Reynolds respectivement. Les resultats sont donnes 
pour une seule pression d’entree (2,12 bars). Us devront 
etre consolides ensuite avec des pressions de travail plus 
elevees. 


4.2. Caracterisation thermique 
du regenerates 

Dans ce travail de mise au point de la methode de 
caracterisation, nous avons effectue les essais de carac¬ 
terisation thermique du mini-regenerateur sous la meme 
pression d’entree que precedemment. Les releves tem- 
porels sont effectues avec une periode d’echantillonnage 
egale a 0,5 s sur une duree de 500 s pour laquelle le 
regime permanent est quasiment atteint aux valeurs du 
nombre de Reynolds considerees. Pour la precision et 
la reproductibilite des experiences, les etalonnages sont 
controles en deux points a l’aide du capteur de reference 
AD22103 : avant de declencher F echelon et a l’equilibre 
thermique. Dans nos conditions operatoires, la calibra¬ 
tion realisee in situ determine une temperature absolue du 
gaz mesuree a 0,3 K pres; soit une incertitude maximale 
inferieure a un pour cent. Sous forme reduite, la tempe¬ 
rature reduite du gaz 0 g est obtenue avec une incertitude 
inferieure a cinq pour cent. En regime permanent, les me¬ 
sures bruitees delivrees par les microthermocouples reve- 
lent un ecart type a egal a 0,1 K. 

La figure 3 illustre les reponses temporelles experi- 
mentales et recalculee par le modele lors de la reponse 
a un pseudo-echelon de temperature du gaz en entree. 
Le profil donne par le modele s’ecarte de celui mesure 
essentiellement aux courts instants puisqu’il prevoit as- 
sez mal l’effet d’accumulation de chaleur des capacites 
thermiques de la matrice et du tube. II s’agit la probable- 
ment d’un probleme lie au traitement numerique des don- 
nees, en particulier a F inversion numerique de la trans- 
formee de Laplace par l’algorithme de Stehfest [8]. Pour 
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2 4 6 8 10 

Temps reduit z = t/R c 


Figure 3. Evolution temporelle des temperatures d’entree et de 
sortie du gaz pour un Reynolds Ret, — 872. 

Figure 3. Time evolution of the entrance and exit gas temper¬ 
atures with a Reynolds number of Ret, = 872. 



Figure 4. Evolution temporelle de la temperature de sortie du 
gaz pour differents nombres de Reynolds. 

Figure 4. Time evolution of the exit gas temperature for 
different Reynolds numbers. 


permettre des comparaisons significatives entre les dif¬ 
ferents releves, nous avons represente sur \& figure 4 les 
courbes de temperatures du gaz en sortie sous forme adi- 
mensionnelle. La normalisation s’effectue au moyen du 
temps caracteristique R c pour la variable temps et des 
temperatures initiale et d’entree du gaz. Au vue des resul- 
tats, on remarque que la configuration generate de l’ecou- 
lement ne subit pas de changement notable puisque l’adi- 
mensionnalisation determine un releve temporel caracte¬ 
ristique de F ensemble des debits concernes. Dans le ta¬ 
bleau III , nous presentons les premiers resultats d’estima¬ 
tion des coefficients d’echange obtenus a partir des me- 
sures de temperatures realisees pour trois configurations 
hydrodynamiques de l’ecoulement (Re h = 462, Re h = 
872 et Re h = 1 060). Le temps caracteristique R c est egal 
a 113,8 s pour Re h = 462, a 61,3 s pour Re h = 872 et a 
51,2 s pour Re h = 1 060. A partir de ces premiers resul¬ 
tats, nous avons constate que le coefficient d’echange h e 
relatif aux pertes externes ne peut pas etre estime en rai¬ 


son de la faible sensibilite des reponses en temperatures a 
ce parametre. En effet, pour ces gammes de temperatures 
du gaz (295-325 K), les echanges avec le milieu am- 
biant sont tres inferieurs aux echanges thermiques entre 
le gaz et la matrice. Les memes difficultes se presentent 
pour l’evaluation du coefficient h ,; en effet, les evolu¬ 
tions temporelles des temperatures sont beaucoup moins 
sensibles aux variations de ce parametre qu’a celles du 
coefficient h s . De plus, les coefficients h e et h t sont etroi- 
tement correles. Les estimations du coefficients h s pre¬ 
sentent des valeurs du meme ordre de grandeur que celles 
donnees par les correlations deja publiees [9-12]. Notons 
que ces correlations existantes ne sont peut-etre pas les 
mieux adaptees a notre configuration et que leur domaine 
de validite respectif n’est pas precise. 

Au vue des resultats, Fhypothese d’une conduction 
thermique axiale negligeable dans la matrice solide par 
les points de contact semble verifiee puisque la valeur 


TABLEAU III 

Comparaison des valeurs estimees et des valeurs publiees des coefficients d’echange. 


R c , s 

Re h 


/i^W-m^-r 1 

h u W-nTAKT 1 

h e , W-m 2 -K 1 

113,8 

Re h = 462 

Valeurs estimees 

20,7 

218 

- 



Correlation [ 10] 

35,7 

- 

10 

61,3 

Ret, = 872 

Valeurs estimees 

42,4 

439 

- 



Correlation [ 10] 

47,6 

- 

9 

51,2 

Re h = 1 060 

Valeurs estimees 

49,2 

448 

- 



Correlation [ 10] 

51,3 

- 

9 
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du nombre de Biot est inferieure a 1 (Bi s ~ 10 -3 ). En 
revanche, cette condition n’est pas verifiee (Bi t >1) pour 
le tube. 

5. CONCLUSIONS ET PERSPECTIVES 

La methode et le montage experimental decrits per- 
mettent de realiser rapidement et avec une precision ac¬ 
ceptable la caracterisation complete d’un mini-regenera- 
teur. Celle-ci passe par Fetablissement des courbes ca- 
racteristiques du frottement et la determination des trois 
coefficients d’echange thermique. L’allure generate de 
la courbe de frottement experimentale correspond assez 
bien a celle obtenue a partir de la relation classique d’Er- 
gun avec cependant des ecarts pour les valeurs extremes 
du nombre de Reynolds (jusqu’a +37,1 %) probablement 
dus a la non-uniformite de la porosite a la paroi de Fen- 
veloppe. Cet aspect devient inevitable et plus manifeste 
dans les mini-systemes. Un effet de la temperature sur 
le frottement, phenomene encore peu etudie, a ete note; 
la perte de charge diminue lorsque la temperature aug- 
mente. 

L’estimation du coefficient externe h e s’est revelee 
impossible, celle de h t assez imprecise en raison de la 
faible sensibilite du modele pour ces parametres dans 
le cadre des conditions experimentales realisees. Pour 
le coefficient /j t , uniquement significatif pour des mini- 
echangeurs, la litterature actuelle ne fournie d’ailleurs 
pas de valeurs de comparaison. L’estimation de /7 S reste, 
quant a elle, coherente avec les relations fournies par les 
auteurs de la speciality, bien que Futilisation de celles- 
ci semble confuse : notations, definitions et domaines de 
validite non precises clairement. 

La. figure 5 represente un prototype de refroidisseur 
a tube a gaz pulse miniature dont la fabrication est is¬ 
sue de la technologie Silicium; sa realisation a ete menee 
par le LPMO (Laboratoire de Physique et Metrologie des 
Oscillateurs) de Besan 5 on. Le regenerates est constitue 
de billes de verre d’une centaine de micrometres de dia- 
metre. L’emploi de la technique decrite doit permettre 
d’optimiser le fonctionnement thermique global de Fap- 
pareil. 
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Figure 5. Refroidisseur a tube a gaz pulse miniature. 
Figure 5. Miniature pulsed tube refrigerator. 
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